Direct observation of a "devil's staircase" in wave-particle interaction 
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We report the experimental observation of a "devil's staircase" in a time dependent system 
considered as a paradigm for the transition to large scale chaos in the universality class of 
hamiltonian systems. A test electron beam is used to observe its non-self-consistent interaction 
with externally excited wave(s) in a Travelling Wave Tube (TWT). A trochoidal energy analyzer 
records the beam energy distribution at the output of the interaction line. An arbitrary waveform 
generator is used to launch a prescribed spectrum of waves along the slow wave structure (a 4 m 
long helix) of the TWT. The resonant velocity domain associated to a single wave is observed, 
as well as the transition to large scale chaos when the resonant domains of two waves and their 
secondary resonances overlap. This transition exhibits a "devil's staircase" behavior for increasing 
excitation amplitude, due to the nonlinear forcing by the second wave on the pendulum-like motion 
of a charged particle in one electrostatic wave. 



PACS numbers: 05.45.Df, 52.40.Mj, 52.20.Dq, 84.40.Fe 



Numerous nonlinear driven systems in physics, 
astronomy, and engineering exhibit striking re- 
sponses with complex phase-locking plateaus 
characterized by devil's staircases. Whereas this 
structure has fascinated theoreticians since its 
discovery, and was observed in various condensed 
matter systems, no dynamical system had so far 
allowed its direct observation. We realized such 
an observation for the motion of a charged par- 
ticle in the potential of two waves. This system 
has already been used as a theoretical and ex- 
perimental paradigm for analysing the transition 
to large scale chaos in conservative systems. Our 
method consists in measuring the velocity distri- 
bution function of a test beam at the outlet of 
a harmonically excited travelling wave tube. A 
devil's staircase is directly measured by varying 
the amplitude of the excitation in perfect agree- 
ment with numerical simulations. This clear ex- 
perimental realization for a theoretical paradigm 
opens the possibility to achieve the detailed ex- 
ploration of chaotic behaviour in non-dissipative 
systems. 



I. INTRODUCTION 

Wave-particle interaction is central to the operation 
of electron devices and to the understanding of basic 
plasma behavior and more generally of systems with long 
range interaction or global coupling. Specifically, travel- 
ling wave tubes are ideal to study wave-particle intcr- 
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action |l|, |^|^. In the experiments reported here, the 
electron beam density is so low that it does not induce 
any significant growth of the waves. When waves are 
launched in the tube with an external antenna, the beam 
electrons can be considered as test particles in the po- 
tential of electrostatic waves. This recently allowed the 
direct exploration of nonlinear particle synchronization 
4] by a single non resonant wave with a phase velocity 
different from the beam velocity. 

In this paper we focus on the resonant case when the 
beam and the waves propagate with similar velocities. 
For a single wave, we observe beam trapping in the poten- 
tial troughs of the wave. When at least two waves propa- 
gate we observe a large velocity spread of the essentially 
mono-kinetic initial beam. This spread can be related 
to a transition to large scale chaos through resonance 
overlap and destruction of invariant KAM (Kolmogorov- 
Arnold-Moser) tori Q ■ We also observe that the velocity 
spreading is not smooth and occurs through successive 
steps as the waves amplitudes increase. 

The motion of a charged particle in two electrostatic 
waves has been used as a paradigm system to describe the 
transition to large-scale chaos in Hamiltonian dynamics 
[y> Q> IE 1^ 113 ■ It is well known that Hamiltonian phase 
space exhibits a scale invariance with infinitely nested 
resonant islands [111. ^^ show that the observed steps 
in the measured velocity distribution function are related 
to the underlying "devil's staircase" behavior. 

Numerous nonlinear driven systems in physics, as- 
tronomy, and engineering exhibit striking responses with 
complex phase-locking plateaus characterized by devil's 
staircase, Arnold tongues, and Farey trees p^ll3lll4Lll5| . 
These structures are mainly seen at work in stationary 
and space-extended objects, when an elastic lattice in- 
teracts with a periodic substrate found in a wide vari- 
ety of condensed matter systems including charge-density 
waves [ig, Josephson junction arrays [i3|, Frenkel- 
Kontorova-type models of friction [l3|, unusual growth 



mode of cr yst als pj| | , and in semiconductor-metal phase 
transition |2(1. They are also at work in the quantum 
Hall effects [2lj . But very few dynamical systems have so 
far allowed their observation. Only indirect observations 
of statistically locked-in transport through periodic po- 
tential landscapes have been made |23|. Here we present 
a direct experimental evidence of a "devil's staircase" for 
wave-particle interaction [23 ■ 

The paper is structured as follows: In Sec. H we present 
the principle of the test beam experiment in the TWT. 
In Sec. Ill, we report the main results on the interaction 
of the test beam with one or two waves. In Sec. IV, we 
briefly recall the main results on the transition to large 
scale chaos in Hamiltonian systems and show some nu- 
merical results for the paradigm motion of a charged par- 
ticle in two waves. In Sec. V, we compare the experimen- 
tal results with the numerical predictions. Conclusions 
and perspectives are drawn in Sec. VI. 



II. DESCRIPTION OF THE TWT 



A travelling wave tube [2J, |25j consists of two main el- 
ements: an electron gun, and a slow wave structure along 
which waves propagate with a phase velocity much lower 
than the velocity of light, which may be approximately 
equal to the electron beam velocity. 

The long travelling wave tube sketched in Fig. 1 has 
been designed to mimic beam-plasma interaction jHQ,!^. 
Its slow wave structure is formed by a helix with axially 
movable antennas. The electron gun produces a quasi- 
mono-kinetic electron beam which propagates along the 
axis of the slow wave structure and is confined by a strong 
axial magnetic field of 0.05 T. Its central part consists of 
the grid-cathode subassembly of a ceramic microwave tri- 
ode and the anode is replaced by a copper plate with an 
on-axis hole whose aperture defines the beam diameter 
equal to 1 mm. Beam currents, /b < 1 /iA, and maximal 
cathode voltages, |K;| < 200 V, can be set independently. 
Two correction coils (not shown in Fig. 1) provide per- 
pendicular magnetic fields to control the tilt of the elec- 
tron beam with respect to the axis of the helix. 

The slow wave structure consists of a wire helix that is 
rigidly held together by three threaded alumina rods and 
is enclosed by a glass vacuum tube. The pump pres- 
sure at the ion pumps on both ends of the device is 
2 X 10~^ Torr. The 4 m long helix is made of a 0.3 mm di- 
ameter Be-Cu wire; its radius is equal to 11.3 mm and its 
pitch to 0.8 mm. A resistive rf termination at each end 
of the helix reduces reflections. The maximal voltage 
standing wave ratio is 1.2 due to residual end reflections 
and irregularities of the helix. The glass vacuum jacket 
is enclosed by an axially slotted 57.5 mm radius cylinder 
that defines the rf ground. Inside this cylinder but out- 
side the vacuum jacket are four axially movable antennas 
labelled from jj 1 to jj 4 and azimuthally separated by 90 "; 
they are capacitively coupled to the helix, each of them 
being characterized by its coupling coefficient which gives 



the ratio of the amplitude of the wave actually excited on 
the helix to the amplitude of the signal launched on the 
antenna; these antennas can excite or detect helix modes 
in the frequency range from 5 to 95 MHz. Only the helix 
modes are launched, since empty waveguide modes can 
only propagate above 2 GHz. These modes have electric 
field components along the helix axis T| . Launched elec- 
tromagnetic waves travel along the helix at the speed of 
light; their phase velocities, v^j, along the axis of the 
helix are smaller by approximately the tangent of the 
pitch angle, giving 2.8 x lO^m/s < v^pj < 5.3 x lO^m/s. 
A beam intensity /b lower than 150 nA ensures that the 
electron beam with mean velocity Vh does not induces any 
significant wave growth over the tube length; since the 
potential drop over the beam radius SV = /b/(47reot'b) is 
negligible, the beam electrons can be considered as test 
particles. 

Finally the cumulative changes of the electron beam 
distribution are measured with a trochoidal velocity an- 
alyzer at the end of the interaction region. A small frac- 
tion (0.5%) of the electrons passes through a hole in the 
center of the front collector, and is slowed down by three 
retarding electrodes. By operating a selection of elec- 
trons through the use of the drift velocity caused by an 
electric field perpendicular to the magnetic field, the di- 
rect measurement of the current collected behind a tiny 
off-axis hole giyes the time averaged beam axial energy 
distribution |2g . Retarding potential and measured cur- 
rent are computer controlled, allowing an easy acquisi- 
tion and treatment with an energy resolution lower than 
0.5 eV. 



III. TEST PARTICLE EXPERIMENTS 

We first checked that, in the absence of any applied 
external signal on the antenna, the test electron beam 
propagates without perturbation along the axis of the 
device. 



A. Trapping of the test beam in a resonant wave 

In the first experiment, we apply an oscillating sig- 
nal at a frequency of / = 30 MHz on antenna ftl. 
According to the helix dispersion relation, a travelling 
wave propagates along the helix with a phase velocity 
v^ — 4.07 X 10^ m/s. Fig. [2t shows a 2D contour plot 
(with logarithmically scaled color coding) of the veloc- 
ity distribution function for a test beam, with intensity 

,„, measured at 



/b = 120 nA and initial velocity fb = 
the outlet of the tube after its interaction, over a length 
of 2.6 m, with the helix mode. It is the result of linear 
interpolation of the recorded measurements obtained for 
different applied signal amplitudes varying from mV to 
1400 mV by steps of 100 mV. Each velocity distribution 
function measurement is obtained by scanning the re- 
tarding voltage of the trochoidal velocity analyzer with a 



step of 61 mV. When the wave ampUtude is gradually in- 
creased, we observe that the width of the velocity domain 
in which the test beam electrons are spread increases 
like the square root of the wave amplitude as expected if 
the electrons are trapped in the potential troughs of the 
wave. No attention must be given to the apparent reg- 
ular discontinuous steps due to bad sampling of the am- 
plitudes and interpolation between successive measure- 
ments. From the measured antenna coupling coefficients 
|27| . the helix wave amplitude (j) can be estimated. The 
resonant or trapping domain in velocity can thus be de- 
duced and is indicated by the continuous lines in Fig. [2^ 
which correspond to w,^ ± 2^/r](f>] here 77 is the charge 
to mass ratio of the electron. We observe a very good 
agreement with measurement. 

Fig.|2t> is obtained in the same way as Fig. [5^ but cor- 
responds to a larger interaction length since the emitting 
antenna is located at L = 3.5 m from the device output. 
Each velocity distribution function measurement is ob- 
tained by scanning the retarding voltage of the trochoidal 
velocity analyzer with a step of 240 mV to speed up data 
acquisition without deteriorating significantly the accu- 
racy of our observations. We observe that the shape of 
the velocity domain in which the test beam electrons are 
spread is very different from Fig. |2t- A new feature ap- 
pearing in Fig. \^ is a further velocity bunching of the 
electrons around their initial velocity for an amplitude 
equal to 600 mV. This phenomenon is also related to the 
trapping of the electrons in the wave. If we refer to the 
rotating bar model _28j to describe the trapped electrons 
motion, we expect oscillations between small spread in 
velocity and large spread in position (corresponding to 
initial conditions for a cold test beam) to large spread in 
velocity and small spread in position after half a bounce 
or trapping period equal to Tb/2 = v^/{2fy/T](j)). To 
the bounce period, we can associate a bounce length 
Lb = v^Tf,. For (j) = 35 mV, we get Lb/2 = 350 cm which 
is precisely the interaction length L of Fig. \Bp. We thus 
confirm that Fig.|2]is displaying the trapping of the test 
beam in a single wave. Varying /b from 10 to 150 nA 
would not change the result since no self-consistent ef- 
fect occurs. One may remark that Fig. [3 also displays a 
systematic accumulation of particles in the high velocity 
region for which wc do not have a simple explanation and 
which probably requires further investigation out of the 
scope of the present paper. 



B. Test beam interacting with two waves 

In the second experiment, we apply, on antenna (j 2 cho- 
sen for its strongest coupling with the helix and located 
at L = 3.5 m from the device output, a signal generated 
by an arbitrary waveform generator made of two compo- 
nents with well defined phases: again one at the same 
frequency of 30 MHz, and a new one at a frequency of 
60 MHz. According to the helix dispersion relation, two 
travelling waves propagate along the helix, the former 



with a phase velocity v^i = 4.07 • 10^ m/s as before, the 
latter with a different phase velocity v^2 = 3.08 • 10^ m/s 
(the fact that the second wave is a harmonic of the first 
one in the laboratory frame is irrelevant in the beam 
frame) . 

Fig. |3Jl (resp. Fig. ^p) is obtained in the same way as 
Fig. [21 for a test beam, with intensity /b = 10 nA and 
initial velocity Vh = v^pi (resp. Vh — Wcp2), whose ve- 
locity distribution function is measured at the outlet of 
the tube after its interaction with the two propagating 
helix modes. As in Fig. [21 the continuous (resp. dashed) 
parabola shows the trapping velocity domain associated 
to the helix mode at 60 MHz (resp. 30 MHz). As it 
appears clearly, the amplitudes of the two modes have 
been appropriately chosen such that these trapping do- 
mains approximately have the same velocity extension 
when taking into account the coupling coefhcient of the 
antenna at the two working frequencies. Both in Fig. [3^ 
and Fig. [3), we observe that, for a threshold ampli- 
tude exceeding the amplitude corresponding to overlap 
of the trapping domains, the velocity distribution func- 
tion exhibits a strong velocity spread. Electrons ini- 
tially trapped inside the potential well of one of the helix 
modes can escape the trapping domain and explore a 
much wider velocity region in phase space. According to 
their initial velocity, electrons can be strongly accelerated 
(Fig. [3^) or decelerated (Fig. ^jp) sticking mainly to the 
border of the new velocity domain they are allowed to 
explore. We shall later relate this behavior to the tran- 
sition to large scale chaos for the motion of a charged 
particle in two electrostatic waves. Such a description 
would suggest that Fig. [^d could be obtained by a mere 
reflection of Fig. [3^ with respect to the average velocity 
of the two waves. The observed asymmetry can result 
from the different transit times of the test particles over 
the finite length of the tube in both cases, and from the 
already noticed accumulation of particles toward higher 
velocities of Fig. [21 



C. Test beam and excitation at a single frequency 

In the third experiment, we apply an oscillating signal 
at a frequency of 30 MHz on antenna ft 2 as in the first 
experiment for maximum interaction length L — 3.5 m. 
But we now consider a test beam, with intensity /b — 
10 nA and initial velocity vi, = 2.7 • 10^ m/s much lower 
than v^, the phase velocity of the helix mode at 30 MHz. 
Fig. 21 is obtained in the same way as Fig. [21 The contin- 
uous parabola indicates the velocity domain associated 
to particle trapping in the potential troughs of the he- 
lix mode at 30 MHz. Since the initial beam velocity 
lies far out of this domain we expect that, for moder- 
ate wave amplitude, the beam electrons will experience 
a mere velocity modulation around their initial velocity, 
with a modulation amplitude increasing linearly with the 
applied signal amplitude. This behavior has been studied 
in detail in Ref. [3] and would generate two main peaks 



around the (oblique straight) continuous hues originating 
in Vh. 

In fact we observe that the electrons spread over a ve- 
locity domain with typical width increasing as the square 
root of the applied signal amplitude as shown by the 
dashed parabola. This strongly recalls the results of 
Fig. [3 Indeed the applied signal generates two waves: 
a helix mode with a phase velocity v^ , and a beam mode 
with a phase velocity equal to the beam velocity Vh ■ The 
beam mode with phase velocity Wb is actually the su- 
perposition of two indistinguishable modes with pulsa- 
tion uj — kvh ± Wb corresponding to the beam plasma 
mode with pulsation Wb = [nhe'^/i'meo)]^''^ for a beam 
with density nb, Doppler-shifted by the beam velocity Vh, 
merging in a single mode since ujh <^ lu in our conditions 
|24|. Thus Fig. 21 shows the test electrons trapping into 
the beam mode. This is confirmed by a careful analysis of 
Fig. ^ which exhibits the same velocity bunching of the 
electrons around their initial velocity as in Fig. |21i, for 
amplitudes obtained by equating the interaction length 
to a multiple of half the trapping length. One also notices 
that the amplitude of the beam mode is lower than the 
amplitude of the launched helix mode; this explains why 
its influence has been neglected in the previous analysis 
of Fig. 01 where two helix modes are externally excited. 
Since, when we increase the applied signal amplitude in 
Fig. ^ the two trapping domains of the helix and the 
beam mode overlap, we observe the same behavior as in 
Fig. 121 Above a certain applied signal amplitude thresh- 
old the distribution function spreads over a much wider 
velocity domain, and electrons can be strongly acceler- 
ated. 

Another striking feature appears in Fig. 21 which is 
best emphasized in the zoom of Fig. \Ej^. The transition 
to large velocity spread does not occur continuously but 
rather occurs by steps when the applied signal amplitude 
increases. As shown in Fig.EJj, plateaus are formed in the 
measured velocity distribution function for the maximum 
interaction length. A closer look at Fig.Oalso reveals the 
presence of such steps in the case of two independently 
launched waves with equal amplitudes. We will now show 
that this generic phenomenon is related to the intrinsic 
structure of Hamiltonian phase space for non integrable 
systems. 



IV. MOTION OF A TEST CHARGED 
PARTICLE AND HAMILTONIAN CHAOS 

We consider charged test particles moving in two elec- 
trostatic waves. The equation modelling the dynamics in 
this case is 



y^ rih<l>i sm{kiX - ujit + ip{). 



(1) 



The motion of the particles governed by this equation 
is a mixture of regular and chaotic behaviors mainly de- 
pending on the amplitudes of the waves |q,|3, n3| and ex- 
hibits generic features of chaotic systems [l^, [l3i u^ US ■ 
Poincare sections of the dynamics (a stroboscopic plot 
of selected trajectories) are displayed in Fig. |S1 Taking 
advantage of the symmetry x' = —x, two different cases 
are superposed in Fig. O The strength of chaos is mea- 
sured by the dimensionless overlap parameter defined as 
the ratio 



(2) 



\v^l - V^2\ 



where 4>i, ki, oji and (fi are respectively the amplitudes, 
wave numbers, frequencies and phases of the two waves; 
rj is the charge to mass ratio of the particle. 



of the sum of the resonant velocity half-widths of the two 
wave potential wells to their phase velocity difference . 

For intermediate wave amplitudes, nested regular 
structures appear as secondary resonances with a wave 
number knm — nki+mk2, a frequency cunm = 'nLJi-\-muj2^ 
a phase (pnm — n(pi + mip2 and an effective amplitude 
<l>nm ~ </'i </'2 with integer n and m, as shown in the 
left side of Fig.JHt for s = 0.55. The self-similar structure 
of phase space results from the infinitely nested higher 
order resonances [ll| which appear in so-called Arnold 
tongues between secondary resonances as predicted by 
the Poincare-Birkhoff theorem |30| . 

For larger wave amplitudes, a wide connected zone 
of chaotic behavior occurs in between the primary 
resonances due to the destruction of the so-called 
Kolmogorov-Arnold-Moser (KAM) tori acting as barri- 
ers in phase space, as shown in the right side of Fig. [^ 
for s — 0.85. Between these two values of s, invariant 
tori are sequentially destroyed and replaced by cantori. 
Transport in velocity (and thereby in {x,v) space) oc- 
curs through the holes of the cantori, which have a fractal 
structure |3l|. When a torus breaks, the flux through it is 
still zero, as the associated turnstile has vanishing area : 
trajectories leak easily through its holes only for s values 
significantly above its destruction threshold. For the dy- 
namics X — — e(sin(x) -I- 0.16ksmk(x — t)), with k — 5/3, 
which is close to our experimental conditions, the thresh- 
old for large scale chaos (destruction of the most robust 
torus) is near s ~ 0.75, as seen from the Poincare section 
of orbits over long times (FigjHl (b-c)). 

If one considers a beam of initially monokinetic parti- 
cles having a velocity equal to the phase velocity of one 
of the waves, this chaotic zone is associated with a large 
spread of the velocities after some time since the parti- 
cles are moving in the chaotic sea created by the overlap 
of the two resonances [3, u^- In analogy with the ex- 
perimental conditions, we compute (with a second order 
centered symplectic integrator) the trajectories of 1 mil- 
lion particles initially uniformly spread in position with 
the same velocity (equal to 1). As seen in fig. [3 as the 
overlap parameter increases, the beam wiggles more and 
is spread over a wider velocity range JS^I . 

As shown in Figs.JHland 9, the transition to large scale 
chaos occurs by vertical steps, related to the presence of 
the nested secondary resonances, with height scaling as 



the amplitudes (jinm- Therefore the border of the velocity 
domain over which the beam is spread exhibits a devil's 
staircase-like behavior. Fig. |Ht is the result of applying 
the same treatment as applied to experimental data in 
section llTTl to the computed velocity distribution function 
obtained for different overlap parameters varying from 
to 1 by steps of 0.01 keeping the wave amplitude ratio 
constant and approximately equal to the experimental 
value of Fig. 0] Each velocity distribution function is 
obtained by recording the velocity histogram of our 10^ 
particles over 1600 boxes with normalized velocity width 
equal to 0.001 after an integration time r equal to 10 
Poincare periods. Even for an overlap parameter equal 
to 1, during the integration time r, the particles could 
not explore the whole chaotic sea between the two main 
resonances. But the particles could typically experience 
one complete trapping motion as shown by the periodic 
velocity bunching and related same horizontally wedge- 
shaped regions as ex per imentally observed in Fig. [3 An 
asymmetric plateau |33l | is formed in the velocity distri- 
bution function as exemplified in the insert of Fig. |SJi. 
The detailed analysis of the border of the distribution 
function displayed in Fig. ^p shows that this plateau can 
be associated to the secondary resonance characterized 
by integers m ~ n ^ 1. 

Fig.l^is obtained in the same way as Fig.|Slbut corre- 
sponds to a longer integration time r equal to 20 Poincare 
periods. Now the particles can explore the whole chaotic 
sea in between the two main resonances, and finer details 
of the "devil's staircase" associated to the presence of the 
nested secondary resonances can be observed. In Fig. ^ 
showing the low velocity border of the computed chaotic 
domain, the position of some of the most important sec- 
ondary resonances is indicated with a label corresponding 
to their characteristic integers (n, m). It appears clearly 
that the presence of a secondary resonance results in a 
step in the 2D contour plot of the distribution function. 
The most striking feature is the great similarity with the 
above experimental results which can now be reexamined 
in the light of the two-waves model. 



V. DISCUSSION 

To allow easier comparison with numerical simulations, 
we compute the overlap parameter of the helix and the 
beam modes. The wave amplitude can be estimated by 
fitting the trapping domain with a parabola as shown 
in Fig. 121 to Fig. 01 As already mentioned, we checked 
that, for the helix mode, this estimate is in very good 
agreement with the amplitude deduced from measuring 
the antenna coupling coefRcient with three antennas |27| . 
Then we display the upper velocity domain border (i.e. 
the location of the upper steep edge of the plateau as 
shown in Fig. [SJd) as the applied signal amplitude is var- 
ied. We thus obtain Fig. \^. The experimentally ob- 
served steps are clearly related to the presence of sec- 
ondary resonances labelled (l,m) with to = 1, 2, 3. 



For a different initial beam velocity, Fig. llOb shows 
the upper velocity domain border as the applied signal 
amplitude is varied by steps of 1.5 V; only the steps 
associated to the main secondary resonances (l,n) for 
n = 1...6 can be observed. The points plotted are the 
outcome of a single set of measurements; the solid line 
is a guide to the eye. By reducing the signal amplitude 
steps to 0.5 V, more details of the devil's staircase can 
be observed and higher order resonances show up. For 
example, in Fig. llOb . more steps are observed such as 
(2,3) = (1,1) + (1,2), (2, 5) and (2, 7) which stand among 
the next strongest resonances in the Farey tree structure. 
In this latter plot, two new, independent measurements 
have been reported (as circles and squares) [l^ to give 
an estimate of the experimental errors, and the solid line 
interpolates between their averages. The behavior re- 
ported in Fig. llOb and Fig. llOb is typical of self-similar 
phenomena. 

It is worth mentioning that we focused our attention 
on the velocity domain in between the two main reso- 
nances. As shown by the upper part of Fig. 2tj we also 
observe vertical steps associated with resonances (1,-5) 
and (1,-6); this feature is also obvious in the numerical 
simulation of Fig. O where a step associated with reso- 
nance (1,-3) clearly appears on the upper part. Since 
numerical simulations deal with a restricted number of 
particles (typically 1 million), only the steps associated 
with lower to's appear. 



VI. CONCLUSION AND PERSPECTIVES 

We have thus exhibited that measuring the energy dis- 
tribution of the electron beam at the output of a travel- 
ling wave tube allows a direct experimental exploration of 
fractal features of the complex Hamiltonian phase space. 
It constitutes a direct experimental evidence of a "devil's 
staircase" in a time-dependent system. 

This knowledge opens new tracks in the study of these 
systems: 

-It allowed accurate experiments aimed at investigat- 
ing Hamiltonian phase space [3J and testing new methods 
of control of Hamiltonian chaos [33 ■ 

-It is a step in the direction of experimental assess- 
ment of the still open question of quasilinear diffusion in 
a broad spectrum of waves |36j | excited with the arbitrary 
waveform generator p37l | . 

-Time resolved measurements of the evolution of elec- 
tron bunches i nject ed with a prescribed phase with re- 
spect to a wave [33| should allow to explore more details 
about the test particle dynamics. 

-The influence of self-consistent effects could be studied 
by gradually increasing the beam intensity. 

Beside their fundamental importance, these studies 
can find direct useful applications in the fields of : 

-plasma physics, 

-control of complex systems. 



-electron devices such as travelling wave tubes and free 
electron lasers where improvement of performances are of 
crucial importance. 
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FIG. 1: Travelling Wave Tube rendering: (1) helix, (2) elec- 
tron gun, (3) trochoidal analyzer, (4) antenna, (5) glass vac- 
uum tube, (6) slotted rf ground cylinder, (7) magnetic coil. 
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FIG. 2: 2D contour plot of the measured velocity distribution 
function (with logarithmically scaled color coding) of a test 
beam (Jb ~ 120 uA, Vh = 4.07 x 10® m/s) trapped in a single 
wave at 30 MHz with trapping domain (continuous curve) 
for increasing amplitude and fixed interaction length: a,) L — 
2.6 m, b) L = 3.5m. 
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FIG. 3: 2D contour plot of the measured velocity distribution 
function (with logarithmically scaled color coding) of a test 
beam (/b = 10 uA) interacting with two waves at 30 MHz 
and 60 MHz with trapping domains (continuous and dashed 
curves) for increasing amplitude and fixed interaction length 
L = 3.6 m: a) Ub = 3.55 x lO*' m/s, b) «b = 4.07 x lO*' m/s 
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FIG. 4: 2D contour plot of the measured velocity distribution 
function (with logarithmically scaled color coding) of a test 
beam (lb ~ 10 uA, Vh = 2.7 x 10 m/s) with applied signal at 
30 MHz for increasing amplitude and fixed interaction length 
L = 3.6 m (trapping domains of helbc and beam modes are 
indicated by contirmous and dashed parabolas). 
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FIG. 5: a) Zoom of Fig. 0] b) Beam velocity distribution 
function for amplitude corresponding to the dashed vertical 
line in the above figure, c) Normalized upper velocity frontier 
of f{v) versus overlap parameter s; velocity normalization is 
such that V — (resp. 1) stands for the helix (resp. beam) 
mode phase velocity; secondary resonances (n, m) are indi- 
cated at velocities mii/{n + van) with k = v^/vh- 
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FIG. 6: Poincare surface of section for the dynamics given by 
X = — £(sin(a;) + 0.16fcsinfc(a; — i)), with k — 5/3, for which 
s = 2.8V£- (a) Left half for e ~ 0.39 (s = 0.55) exhibits 
island chains of secondary resonances at rational velocities 
m/{n + m), as seen from 15 orbits which do not mix ; right 
half for e ~ 0.92 (s — 0.85) exhibits large scale chaos as seen 
for two orbits, (b) Two orbits related to upper and lower 
main chaotic domains for s = 0.735, each iterated over 10'' 
Poincare periods ; to allow an easy comparison, only points 
with — 37r < (x mod Gtt) < are plotted for one orbit, and 
< (a; mod Gtt) < 37r for the other orbit, (c) Similar plot for 
s = 0.74. 
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FIG. 7: Evolution of the beam initially at w = 1 for the same 
dynamics as Fig. El for values of s around the most prominent 
step at times r = 67r/5 (bold, red dots) and r = 12-k (small, 
black dots). 
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FIG. 8: Numerical devil's staircase, a) Probability distribu- 
tion function f(v) (with logarithmically scaled color coding) 
of the velocity for the dynamics of Fig. 15] as a function of 
overlap parameter s = 2.?>^/e after a time r = 12t7, i.e. 10 
Poincare stroboscopic periods (inset, f{v) for the value of s 
indicated by the vertical dashed line), b) Velocity frontiers of 
f{v) versus s for t — 127r; dashed oblique lines starting from 
V — and v — 1 indicate the primary resonances trapping do- 
mains; secondary resonances (n, m) are indicated at rational 
velocities mk/{n + mk). 
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FIG. 9: Numerical devil's staircase after a time r = 247r, i.e. 
20 Poincare stroboscopic periods, a) Probability distribution 
function f{v). b) Velocity frontiers of f{v). 



16 



0.95 




0.2 0.4 0.6 O.S 

overlap parameter 



0,95 




0.65 



0,4 0.6 O.S 

overlap parameter 



FIG. 10: Experimental devil's staircase: plot versus overlap 
parameter of the upper normalized velocity frontier of /(«) 
for a test beam (7b = 10 uA, Vh ~ 2.5 x 10^ m/s) with applied 
signal at 30 MHz for increasing amplitude and fixed inter- 
action length L = 3.6 m with (a) coarse and (b) fine signal 
amplitude scan; velocity normalization is such that v — 
(resp. 1) stands for the helix (resp. beam) mode phase ve- 
locity; secondary resonances (n, m) are indicated at velocities 
mn/{n + mn) with k = Vi^/vb. 



